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Abstract — Rainfall-Runoff simulation of any basin plays an 
important role for any civil engineering project and proper 
planning and management of the water resources. Runoff 
estimation in catchment is a serious challenge for the hydrologist 
where most of the watersheds are ungauged. Various 
approaches for runoff estimation are available ranging from 
lumped to physically based distributed models. This paper 
describes the use of RRL AWBM (Rainfall Runoff Library 
Australian Water Balance Method Model), to investigate its 
performance, efficiency and suitability in Shipra river basin in 
Madhya Pradesh, India. The AWBM is a catchment water 
balance model capable of simulating runoff using daily or hourly 
rainfall data. AWBM model was developed using daily weighted 
rainfall, evapotranspiration and observed runoff time series of 
11 years period from 1996 to 2006. The model was calibrated for 
the six years period from 1996 to 2001 and validated for five 
years period from 2002 to 2006. The reliability and 
performance of AWBM model was evaluated based on Accuracy 
criteria such as Coefficient of determination (R 2 ), Nash-Sutcliffe 
Efficiency Index (El), Root Mean Square Error (RMSE), and 
correlation coefficient (r). The coefficient of determination R 2 
value for calibration and validation period was observed 0.656 
and 0.496 respectively indicating good agreement between the 
observed and simulated runoff. The Nash-Sutcliffe Efficiency 
(El) for calibration and validation is 65.40% and 48.40% 
respectively. The model was found suitable for Shipra basin in 
predicting daily runoff with good degree of accuracy. 

Index Terms — AWBM, RRL, Rainfall runoff modeling, 
Shipra basin 

I. Introduction 

Water is one of the most important natural resources 
available to mankind, which is primarily responsible for 
sustenance of life on the earth. Rainfall-Runoff is a major 
component of the hydrologic cycle plays a key role in water 
resources management. A rainfall-runoff model is a 
mathematical representation describing the rainfall -runoff 
relations of a catchment area, drainage basin or watershed. 
Rainfall-runoff models are classified as Deterministic, 
Stochastic, Conceptual, Theoretical, Black box, Continuous, 
Event, Complete, Routing and Simplified [1] (Linsley, 1982). 
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Because of its non-linear and multi-dimensional nature, 
rainfall-runoff modeling is extremely complicated [2] 
(Lipiwattanakarn et al., 2004). Many hydrologic models are 
available; varying in nature, complexity and purpose [3] 
(Shoemaker et al., 1997). The widely known rainfall -runoff 
models identified are the rational method, Soil Conservation 
Services (SCS) Curve Number method, and Green-Ampt 
method [4] (Galkate et al., 2011). 

The AWBM stands for Australian Water Balance Model is a 
part of Rainfall-runoff library (RRL) developed by 
Cooperative Research Centre for Catchment Hydrology 
(CRCCH), Australia [5]. It is a computer based conceptual 
rainfall runoff model. Number of studies has been conducted 
using RRL AWBM model. Sharifi, 1994 [6] developed 
AWBM model for the small catchment of New England in 
Australia. The AWBM model was used for rainfall runoff 
modelling and found suitable in predicting daily runoff. 
Boughton, 2006 [7] worked on AWBM model for the 
hydrologic rainfall runoff modelling and concluded that this 
model gave more than two-third good calibrated values and 
scope of improvement of the data sets is present for further 
study. Kumar, 2013 [8] simulated runoff of Tadepalli mandal, 
Guntur, Andhra Pradesh using RRL toolkit and obtained good 
co relation of 0.76 between observed and simulated runoff. 
Balvanshi and Tiwari, 2015 [9] developed a Rainfall Runoff 
model using AWBM for runoff simulation in Bina river basin 
of Madhya pradesh. Jamal, 2011 [10], Haque, 2014 [11] and 
Yu, 2015 [12] also used AWBM model for rainfall runoff 
simulation. The main objective of this study was to develop 
rainfall runoff model for runoff simulation using RRL toolkit 
AWBM model for Shipra river basin of Madhya Pradesh in 
India. 

II. STUDY AREA 

Shipra is one of the sacred rivers in Hinduism. The holy city 
ofUjjainis situated on its right bank. After every 12 years, 
the Kumbh Mela (also called Simhastha) takes place at Ujjain 
on the city's elaborate riverside Ghats, besides yearly 
celebrations of the river goddess Kshipra. Shipra river basin 
has been extended between 76° 06 ' 20 " and 75° 55 ' 60 " North 
Latitude and 22° 97' 00' 'and 23° 76' 20' ' East Longitude and 
covers area of 5679 sq. km. The river traverses total course of 
about 190 km. It originates from Kakribardi hills of Dhar and 
flows north across the Malwa Plateau to join the Chambal 
River. The average annual rainfall of area is about 931.87 
mm. The rainfall in the area is due to the southwest monsoon 
which starts from the middle of June and ends in last week of 
September. The topography is generally rolling to undulating. 
The index map of Shipra river basin is shown in Fig. 1 . 
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III. METHODOLOGY 

A. About AWBM RRL model 

The Cooperative Research Centre for Catchment 
Hydrology (CRCCH), Australia developed the 
Rainfall-runoff library (RRL) AWBM model. It is available 
in public domain. AWBM stands for Australian Water 
Balance Model. This model is a conceptual, lumped 
rainfall-runoff model. AWBM is a catchment water balance 
model that develops relationship between rainfall and runoff. 
The AWBM takes daily time series rainfall and 
evapotranspiration data as input and the output we get is 
daily catchment runoff. Several calibration optimisers and 
display tools are featured in this AWBM model to aid model 
calibration. There are total of nine parameters involved in 
the model structure. The model uses 3 surface stores to 
simulate partial areas of runoff. Al, A2 and A3 are the 
partial area of smallest, middle and largest store 
respectively. Cl, C2 and C3 are the capacity of smallest, 
middle and largest store respectively. The water balance of 
each surface store is calculated independently of the others. 
AWBM estimates moisture balance at either daily or hourly 
time steps. The rainfall is added to each of the three surface 
stores and simultaneous subtraction of evapotranspiration is 
done from each stores. If evapotranspiration is greater than 
the moisture content then value becomes negative and the 
value is set to zero. If moisture in the store becomes greater 
than its capacity then runoff is generated. When runoff 
occurs from any store, part of the runoff becomes recharge of 
the base flow store. The fraction of the runoff used to 
recharge the base flow store is BLI runoff. BLI is the base 
flow index which is the ratio of base flow to total flow in the 
stream flow. The remains of the runoff, i.e. 
(1.0-BLI)*runoff, becomes the surface runoff. The base flow 
store runs down at the rate of (1.0 - K)*BS, where BS is the 
current moisture in the base flow store and K is the base flow 


recession constant. The surface store operates in the similar 
way as the base flow store, and runs down at the rate of 
(1.0-KS)*SS, where SS is the current moisture in the surface 
runoff store and KS is the surface runoff recession constant. 
The schematic diagram of the AWBM model is shown in 
Fig. 2. 



CI-C3 = Surface Storage Capacities Jiiwncc 

AI-A3 = Partial Areas Represented by Surface Storages 
BFI = Bnscflo\\ Index 
K = Daily Bascflow Recession Constant 
BS = Current Volume in Bascllow Store (*) 

KS = Daily Surface Flow Recession Constant 
SS = Current Volume in Surface Routing Store (*) 

Lig. 2: Structure of the AWBM model 

B. MODEL SETUP 

The RRL AWBM model was setup to carry out 
rainfall-runoff modeling in Shipra river basin at Ujjain G/d 
site having catchment area 2102 km 2 and average annual 
rainfall 931.87 mm. The input information of daily rainfall, 
runoff and evapotranspiration for the period of 1 1 years from 
1996 to 2006 was converted to tts format and used for model 
development. 

C. MODEL CALIBRATION 

Calibration is a process to standardize estimated or 
simulated values by using deviations from observed values 
for a particular basin. It thus helps in deriving correction 
factors that can be applied to generate predicted values. 
These simulated values are consistent with the observed 
values. When the AWBM model was set up, model was 
calibrated from 1st Jan 1996 to 31st Dec 2001. The model 
was first run in auto-calibration mode using genetic 
algorithm. The model parameters were optimized manually 
to obtain best set of model parameter simulating runoff with 
high accuracy. 

D. MODEL VALIDATION 

Model validation means forjudging the performance of the 
calibrated model over the portion of historical records which 
have not been used for the calibration. RRL AWBM model 
thus calibrated was then validated for the remaining period 
of five years from 2002 to 2006. During validation the set of 
model parameters obtained during the calibration was used. 
The statistics of the simulated results were analyzed and 
outputs of the model were checked to compare the simulated 
and observed runoff to verify the capability of calibrated 
model to simulate the runoff. 
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E. ACCURACY CRITERIA 

Accuracy of the model can be examined on the basis of 
Coefficient of Determination (R 2 ), Efficiency Index (El) and 
Root Mean Square of Error (RMSE). The use of the 
coefficient of determination (R 2 ) is to test the goodness of fit 
of the model and explains predictive power of model. 

r2 = SILifap - VoXv* ~ 

- 5 *) 2 ] csr =1 (<7, - ^> 2 ] 

Where, q 0 = observed flow, q 0 = mean value of observed flow, 
q s = simulated flow and n = number of data points. 

The reliability of the model was evaluated on the basis of 
Efficiency Index (El) as described by the Nash and Sutcliffe 
[13] 



Where, q 0 = observed flow, q 0 = mean value of observed flow, 
q s = simulated flow and n = number of data points. 

While analyzing the accuracy of the model, the objective 
function was to minimize the Root Mean Square of Error 
(RMSE) between the observed and simulated runoff. 

RMSE = ^Z” =1 (q c -q s ) 2 

Where, q 0 = observed flow, q 0 = mean value of observed flow, 
q s = simulated flow 

F. SENSITIVITY ANALYSIS 

The sensitivity analysis is done so as to understand how 
sensitive a model is to certain parameters. This is useful to 
identify with how the model functions and also what 
parameters need more consideration than others. If the model 
is extensively affected by a particular parameter than the 
focus of calibration should be on that parameter. The AWBM 
provides a feature to examine the sensitivity of all the model 
parameters. 


IV. RESULTS AND DISCUSSION 

The RRL AWBM model was setup to carry out 
Rainfall-runoff modeling in Shipra river basin at Ujjain G/d 
site having catchment area 2102 km 2 . The AWBM model was 
calibrated for six years period from 1996 to 2001 and then 
validated for the remaining period of five years from 2002 to 
2006. Optimized parameters of RRL AWBM model found 
using genetic algorithm. The graph presenting comparison 
between observed and simulated discharge during model 
calibration is shown in Fig. 3 which gives the idea and view of 
best match obtained during the model calibration. Fig. 4 
shows the comparison between observed and simulated 
discharge during the calibration of AWBM model for specific 
year 1998 indicates the very good match. Similarly Fig. 5 and 
Fig. 6 show the comparison between observed discharge and 
simulated discharge during the validation of AWBM model. 
The observed and simulated hydrographs were found 
matching well for peak and low flows as well reasonably. 


Table 1: Optimized parameters of RRL AWBM model 


Parameter 

Optimized 
Parameter Values 

Range of the 
Parameters 

Ai 

0.134 

0 . 000 - 1.000 

A 2 

0.433 

0 . 000 - 1.000 

BFI 

0.494 

0 . 000 - 1.000 

Cl 

29.85 

0-50 

C 2 

116.90 

0-200 

C 3 

193.67 

0-500 

K-base 

1.00 

0 . 000 - 1.000 

Ksurf 

0.252 

0 . 000 - 1.000 



Fig. 3: Comparison between observed and simulated 
discharge for calibration 



Fig. 4: Comparison between observed discharge and 
simulated discharge during the calibration of AWBM model 
for specific year 1998. 



Fig. 5: Comparison between observed and simulated 
discharge for validation 
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Fig. 6: Comparison between observed discharge and 
simulated discharge during the validation of AWBM model 
for specific year 2003 


The Performance of AWBM model was evaluated based on 
Accuracy criteria such as Coefficient of determination (R 2 ), 
Nash-Sutcliffe Efficiency Index (El), Root Mean Square 
Error (RMSE), and correlation coefficient (r). The coefficient 
of determination R 2 value for calibration and validation 
period is 0.656 and 0.496 respectively which indicating good 
agreement between the observed and simulated runoff. The 
Nash-Sutcliffe Efficiency (E) for calibration and validation is 
65.40% and 48.40% respectively. The value of these accuracy 
parameters indicates that the AWBM model is performing 
well in predicting runoff. 


Table 3: Accuracy Parameter values during calibration, 
validation and total period 


Accuracy criteria 

Calibration 

Validation 

Total 

period 

Coefficient of 
determination (R 2 ) 

0.656 

0.496 

0.622 

Nash Sutcliff 
efficiency (El) % 

65.40 

48.40 

60.40 

Correlation 
coefficient (r) 

0.810 

0.790 

0.789 

Root mean square 
error (RMSE) 

58.40 

56.74 

57.65 


The sensitivity analysis of AWBM model parameters was 
performed so as to understand how sensitive a model is to 
certain parameters. The AWBM provides a feature to 
examine the sensitivity of model parameters. Sensitivity 
analysis indicates that four parameters are sensitive and four 
are non sensitive. 

T abl e 2: Sensitivity analysis of AWBM model para meters 


Model Parameters 

Sensitivity Analysis 

Ai 

Non sensitive 

A 2 

Non sensitive 

BFI 

Sensitive 

Q 

Non sensitive 

c 2 

Non sensitive 

C 3 

Sensitive 

Kbase 

Sensitive 

K sur f 

Sensitive 


V. CONCLUSION 

In this study, AWBM RRL model has been developed and 
tested for the performance and suitability in Shipra river basin 
of Madhya Pradesh, India. AWBM is lumped and conceptual 
rainfall-runoff model developed by CRCCH, Australia and 
available in public domain. AWBM model was calibrated and 
validated using daily weighted precipitation, daily 
evapotranspiration and daily observed runoff time series of 1 1 
years period from 1996 to 2006. The model was calibrated for 
the six years period from 1996 to 2001 and validated for five 
years period from 2002 to 2006. Genetic algorithm was used 
for the optimization of the model parameters of the AWBM 
model. The sensitivity analysis was carried out on all the 
AWBM parameters so as to find the sensitive parameters. The 
AWBM model was tested evaluated based on Accuracy 
criteria such as Coefficient of determination (R 2 ), Efficiency 
Index (El), Root Mean Square Error (RMSE), and correlation 
coefficient (r). The coefficient of determination R 2 value for 
calibration and validation period is 0.656 and 0.496 
respectively which indicating good agreement between the 
observed and simulated runoff. The Nash-Sutcliffe 
Efficiency (E) for calibration and validation is 65.40% and 
48.40% respectively. The model was found suitable for 
Shipra basin in simulating hydrological response of the basin 
to the rainfall and predicting daily runoff with good degree of 
accuracy. 
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